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ABSTRACT 

This  paper  uses  the  opportunity  given  by  the  formalism 
of  geometric  algebras  to  perform  colour  image  processing 
by  describing  one  colour  with  a  geometric  vector  in  the  3D 
space.  The  paper  proposes  to  manipulate  colours  in  the  RGB 
space  where  properties  such  as  saturation  or  luminance  can 
be  described  with  the  help  of  geometric  operations.  These 
geometric  operations  are  applied  on  vectors  encoded  in  ^3. 
A  spatial  filtering  operation  on  images  with  this  model  leads 
to  a  colour  gradient  that  can  be  used  for  example  for  colour 
edge  detection. 

1.  GEOMETRIC  ALGEBRA 

1.1  Introduction 

The  (Si ’n  geometric  algebra  is  the  2n -dimensional  vector  space 
extended  with  the  geometric  product  from  the  yn  which  is  a 
^-dimensional  vector  space  of  real  numbers.  The  geometric 
algebra  manipulates  entities  called  multivectors  that  are  the 
extension  of  vectors  to  higher  dimension.  They  are  graded- 
algebras,  that  is  to  say  scalars  are  represented  in  geometric 
algebra  by  zero-dimensional  quantities  (0-graded)  for  exam¬ 
ple.  Vectors  are  one-dimensional  directed  quantities  (also 
represented  by  arrows),  so  1 -graded.  In  geometric  algebra, 
there  are  grade-2  entities  termed  bivectors  which  are  plane 
segments  (for  example  circles)  endowed  with  orientation.  In 
general  a  ^-dimensional  entity  is  known  as  a  k-  vector.  For  an 
overview  on  geometric  algebras  see  [3,  6,  7]  for  example. 

In  this  paper  we  limit  the  focus  on  ^3  algebra  because  we 
use  colour  as  vectors  in  the  3D  space.  In  this  algebra  the  ba¬ 
sis  (^0,^1,^2^3,^23^31^12,^123)  includes  a  scalar,  vectors, 
bivectors  and  a  trivector  (3-graded)  also  called  pseudo-scalar. 
The  vectors  are  the  same  as  those  in  the  basis  (e  1,^2, £3) 
of  y3.  The  bivectors  ^23,  £31  and  e\2  are  the  geometric 
product  of  two  vectors,  ^2^3,  £3^1  and  e\e2  respectively. 
The  remaining  parts  are  the  scalar  eo  =  1  and  the  trivector 
£123  =  £i£2^3-  Any  multivector  A  =  Ao^o  +A\e\  +  A2^2  + 
A3^3  -\-A23e23  +A3i^3i  +Ai2^i2  +Ai23^i23  G  ^3  is  a  linear 
combination  of  the  basis  entities  also  called  blades.  The  k- 
grade  part  of  a  multivector  A  is  given  from  the  grade  op¬ 
erator  (A)k.  For  instance  here,  the  2-grade  part  of  A  is 

(A)  2  =  A23^23  +A3i<?3i  +Ai2^12- 

Colour  image  processing  can  be  done  using  different 
colour  spaces.  Most  devices  as  screens  for  example  use  the 
RGB  space  to  display  colours.  In  order  to  avoid  numeric 
approximations  inescapable  from  conversions  between  dif¬ 
ferent  colour  spaces,  we  choose  to  work  directly  in  RGB. 
This  colour  space  will  be  considered  with  an  Euclidian  met¬ 
ric  to  allow  geometric  calculations.  CIE  L*a*b*  colour  space 
respects  the  human  characteristics  to  distinct  colours  with 


Euclidian  metric  so  to  perform  Euclidian  calculations  from 
RGB  is  an  approximation.  We  propose  to  encode  a  N  x  M 
colour  image  /  by  its  three  colour  components  on  the  three 
vector  parts  of  a  multivector  in  %  :  J^Tq1  T^=o  /1  (xA)e\  + 
f2(x,y)e2+  /3(x,y)e3.  f\{x,y),  f2(x,y)  and  /3 (x,y)  are  the 
red,  green  and  blue  components  of  the  pixel  at  coordinates 

(x,y)- 

1.2  Products 

Products  in  geometric  algebras  are  the  elementary  opera¬ 
tions: 

•  Outer  Product:  This  product  is  also  called  exterior  pro¬ 
duct,  denoted  by  A.  It  generates  the  subspaces  defined  by 
linearly  independent  combinations  of  blades.  For  exam¬ 
ple,  multiplying  two  linearly  independent  vectors  gives  a 
bivector.  If  this  bivector  is  again  multiplied  by  a  linearly 
independent  vector,  the  result  is  the  pseudo- scalar. 

•  Inner  Product:  This  product  also  called  interior  product 
or  left  contraction,  denoted  by  J ,  is  used  to  give  the  notion 
of  orthogonality  between  two  multivectors  [2] .  Let  A  be 
a  a- vector  and  B  be  a  ^-vector,  then  AJB  is  a  (b  —  a) - 
vector  subspace  of  B  and  orthogonal  to  subspace  A.  If 
b  <  a  then  AJ  B  =  0.  Two  vectors  are  also  orthogonal  if 
their  inner  product  is  null. 

•  Geometric  Product:  This  product  is  an  associative  law 
and  distributive  over  the  addition  of  multivectors.  In  ge¬ 
neral,  the  result  of  the  geometric  product  is  a  multivector. 
If  used  on  1 -vectors  a  and  b ,  this  is  the  sum  of  the  inner 
and  outer  product:  ab  =  a\b  +  a  Ab.  Note  also  that  the 
geometric  product  is  not  commutative. 

•  Scalar  Product:  This  product  also  called  dot  product,  de¬ 
noted  by  •,  is  used  to  define  distances  and  modulus. 

1.3  Definitions 

Let  A  be  a  multivector  of  <$3. 

•  The  reversion  of  A,  denoted  A  is  defined  by  A  =  (A)o  + 
(A)i  —  (A) 2  —  (A) 3.  Reversion  is  needed  to  define  Norm 
and  Inverse. 

•  Norm  or  Modulus  of  a  multivector  A  is  defined  by  the 
scalar  product  of  A  by  its  reversion  :  |  A|  =  A  •  A  = 

\f  (AA)o 

Let  now  B  in  & 3  be  the  geometric  product  of  n  1 -vectors 
(n  G  N),  B  is  called  a  versor. 

•  Every  versor  has  an  inverse  which  means  for  every  B, 
there  is  a  B-1  asserting  BB-1  =  1;  the  inverse  B_1  of 

B  equals  to  . 

BB 

•  Note  that  every  multivector  is  not  invertible. 
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2.  GEOMETRIC  TRANSFORMATIONS 

Since  the  colour  information  of  an  image  is  encoded  in  the 
vectorial  part  of  a  ^3  multivector,  we  can  manipulate  a  colour 
in  the  3D  space  with  the  geometric  transformations  given  by 
the  algebra  formalism  (cf.  Figure  7): 

•  Translation:  Let  vi,V2,v?)  E  £f33,  the  vi  translation  by  V2 
is  given  by  their  sum  vt  =  v\  +  V2. 

•  Projection:  Let  (vi ,  V2,  vii)  E  £f3,  the  v\  projection  over  V2 
is  given  by  :  vy  =  (vij  V2)v2~1. 

•  Rejection:  Let  (vi,V2,vx)  E  Sf33,  the  rejection  of  v\  with 

respect  to  V2  is  given  by  :  =  (vi  A  V2)v^ l. 

•  Reflection:  Let  (vi,V2,vr)  E  £f33,  the  reflection  vr  of  v\ 
with  respect  to  V2  is  given  by  :  vr  =  V2V1V2  x. 


V-L 


Figure  1:  Basis  geometric  transformations  from  two  ^3 
vectors  v\  and  V2  :  vt  is  the  v\  translation  by  V2  ;  vn  is  the 
projection  of  v\  on  V2  ;  is  the  rejection  of  v\  with  respect 
to  V2  and  vr  is  the  reflection  of  v\  with  respect  to  V2. 


3.  COLOUR  GRADIENT 

From  the  study  of  the  geometric  transformations  seen  before 
with  the  vectors  of  ^3,  we  propose  to  manipulate  the  colours 
of  images  as  vector  in  the  3D  colour  space.  Felsberg  used  ge¬ 
ometric  algebras  to  characterize  the  luminosity  variations  in 
greyscale  images  by  a  2D  analysis  [5].  A  first  application  of 
geometric  algebras  with  colour  images  was  given  by  Schlem- 
mer  [11].  In  his  paper,  Schlemmer  decomposed  images  in  the 
LUV  colour  space  and  merged  two  gradients.  The  first  one  is 
given  by  a  Canny  method  on  the  luminosity  part  of  the  LUV 
colour  space.  The  second  gradient  is  calculated  with  the  two 
chromatic  parts  U  and  V  with  a  2D  gradient  produced  from 
the  vector  field  theory  using  the  geometric  algebra 

We  propose  in  this  paper  to  use  the  formalism  of  ^3,  that 
is  to  say  to  use  the  3D  information  colour  vectors,  to  perform 
a  colour  edge  detection.  In  fact,  in  the  homogenous  regions 
of  a  natural  colour  image,  pixels  represented  by  vectors  are 
close  to  their  neighbours  in  the  colour  space.  A  colour  edge 
can  then  be  detected  by  a  discontinuity  in  this  neighbour¬ 
hood.  Geometric  manipulations  on  the  colour  vectors  of  i- 
mages  can  help  us  to  detect  such  discontinuities. 

A  first  approach  was  given  by  Sangwine  in  [9].  In  his 
paper,  Sangwine  proposed  to  use  the  formalism  of  quater¬ 
nions  to  encode  the  colour  information  of  an  image.  This 


formalism  is  also  able  to  manipulate  vectors  in  the  3D  space 
but  has  limits.  With  the  definition  of  the  convolution  on  a 
quaternionic  colour  image  in  [8,  10],  he  applied  two  specific 
conjugate  quaternionic  filters  on  an  image.  This  fulfilled  an 
average  vector  of  every  pixel  in  the  filters  neighbourhood  re¬ 
flected  with  respect  to  the  greyscale  axis.  This  reflection  ope¬ 
ration  gives  a  comparison  of  the  colour  vectors  in  the  neigh¬ 
bourhood  of  the  filters  (dimension  of  filters  is  3  x  3).  The 
filtered  image  ( cf.  Figure  2)  is  a  greyscale  image  almost  ev¬ 
erywhere,  because  in  homogeneous  regions  the  vector  sum 
of  one  pixel  to  its  neighbours  reflected  by  the  grey  axis  has 
a  low  saturation.  Saturation  is  the  distance  between  a  colour 
vector  and  the  greyscale  axis.  This  is  the  case  of  <74  +  Jdqsjl 
for  instance  in  Figure  3  where  the  quaternion  /i  represents 
the  greyscale  axis,  Jl  its  conjugate  and  /K73/I  the  reflection  of 
#3  with  respect  to  jLl  [9] .  However,  if  you  compare  pixels  like 
q\  and  q2  for  example  in  Figure  3,  you  should  detect  a  colour 
edge  as  they  are  in  colour  opposition.  Therefore,  they  present 
a  vector  sum  far  from  the  grey  axis.  Edges  are  thus  coloured 
due  to  this  high  distance.  The  limit  of  this  method  is  that 
the  filters  can  be  applied  horizontally  from  left  to  right  and 
from  right  to  left  for  example  but  without  giving  the  same  re¬ 
sults.  In  Figure  3,  the  vector  sum  q2  +  /igi/l  is  different  from 
the  vector  sum  q\  +  \lq2F  produced  when  the  filters  are  ap¬ 
plied  from  opposite  side  (ie.  leftwise  or  rightwise).  An  edge 
with  the  Sangwine  method  will  be  detected  by  full  coloured 
pixels  in  opposition  to  homogenous  regions  which  will  be 
represented  by  grey  pixels. 


(a)  (b) 


Figure  2:  Sangwine  edge  detector  result  :  (a)  original  im¬ 
age;  (b)  zoom  on  the  hat  edge  of  the  (a)  picture  after  applica¬ 
tion  of  Sang  wine’s  filters  from  right  to  left  horizontally. 

In  this  paper,  we  propose  to  detect  edges  by  a  colour  gra¬ 
dient.  The  method  of  Sangwine  gives  us  the  starting  point. 
In  fact,  we  have  a  map  of  every  pixel  reflected  by  the  grey 
axis  and  compared  to  its  neighbours.  The  more  the  vec¬ 
tor  result  is  far  from  the  grey  axis,  the  more  it  is  far  as 
well  from  its  neighbours.  Convolution  product  is  defined 
in  [4]  with  geometric  algebras.  The  filtered  image  is  given 
by  f(x,y)  =  H/(Ti,T2)/((x- Ti),(y- T2))r(Ti,T2)  with 
Sangwine’s  generalized  filters  to  5%: 
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We  suggest  to  determine  the  distance  of  Sangwine  com¬ 
parison  vector  sum  to  the  grey  axis:  its  saturation.  This 
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Figure  3:  Sangwine’s  edge  detector  scheme  :  /u  is  the 

grey  axis;  /igi/I  (resp.  is  the  reflected  vector  of  q\ 

(resp.  ^3)  with  respect  to  /i;  the  comparison  vector  between 
q\  and  q 2  (resp.  between  q 3  and  <74)  is  given  by  q2  +  Liq  1J1 
(resp.  q4  +  {I qsjiy,  q\  +  flqijl  is  near  the  grey  axis  so  the 
colour  seems  grey  but  <72  +  /igi/I  is  far  from  the  grey  axis  so 
Sang  wine’s  filter  has  detected  an  edge  as  this  vector  is  more 
coloured. 


distance  S  is  the  norm  of  the  rejection  of  the  colour  vector 
sum  v  with  respect  to  the  grey  axis  /1.  We  then  get  it  by 
the  formula  of  the  rejection:  S  =  |(v  A/i)/l-1|  (cfi  Figure 
4).  The  same  colours  v\  and  V2  can  give  two  different  colour 
vectors  vsumi  =  V2  +  fivip  and  vsum2  =  vq\  +  flv 2J1  by  the 
Sangwine’s  method.  Note  that  our  approach,  which  gives 
a  saturation  gradient,  is  independent  from  the  path  (leftwise 
or  rightwise)  applied  to  convolute  the  filters  as  the  distance 
Si  =  I (vsumi  A  1  =  | (vSum2  A  This  saturation 

filtering  is  applied  to  the  horizontal,  vertical  and  both  dia¬ 
gonal  directions  by  @3  convolution  product.  The  maximum 
of  these  directional  gradients  is  then  selected  to  make  the  fi¬ 
nal  colour  gradient  filter  by  maximum  saturation.  This  final 
filtering  operation  is  not  linear  only  because  the  ’’maximum” 
operator  interferes. 

As  the  method  described  before  is  a  saturation  gradient,  a 
problem  appears:  discontinuities  in  colours  are  not  just  satu¬ 
ration  ones.  Indeed,  two  different  colours  can  have  the  same 
saturation;  fortunately  they  can  be  disjoined  by  their  differ¬ 
ence  of  luminance. 

To  improve  our  algorithm  we  use  the  opportunity  given 
by  the  embedding  of  the  Geometric  Algebra.  In  fact,  ap¬ 
plying  the  geometric  product:  =  f(x,y)fi  on  every 

pixel  of  the  original  image  gives  us  a  geometric  comparison 
of  the  pixels  to  the  greyscale  axis.  Then  the  result  of  this 
operation  adds  luminance  information.  It  is  the  sum  of  two 
parts:  f(x,y)  A  /I  +/(x,y)J/i.  The  first  part,  f(x,y)  A  fi,  is 
bivectorial  and  allows  to  compare  the  geometry  between  the 
colour  vector  and  /I.  When  this  part  is  null  for  example,  the 
colour  vector  is  collinear  to  the  greyscale  axis.  The  second 
part,/(jc,y)J/i,  is  a  scalar,  and  the  projection  of  the  colour 
vector  over  the  grey  axis  :  its  intensity. 

A  boolean  mask  is  produced  by  thresholding  the  bivec- 


Figure  4:  Saturation  gradient  :  /ivi/I  (resp.  /1V3/I_1)  is 
the  reflected  vector  of  v\  (resp.  V3)  with  respect  to  the  grey 
axis  jU;  our  comparison  vector  between  v\  and  V2  (resp.  be¬ 
tween  V3  and  V4)  gives  the  distance  S  =  |  ( vsum  A  AO H  ‘I  of 
Sangwine’s  vector  sum  vsu„n  =  v2  +  A''’!  A'1  (resP-  vsum2  = 
V4  +  /IV3/T-1)  from  the  grey  axis  (orange  arrows);  An  edge  is 
detected  by  a  high  saturation  like  the  orange  arrow  between 
M  and  + 


torial  part  where  its  norm  is  very  small  to  keep  the  coordi¬ 
nates  where  the  colour  vector  is  mere  collinear  to  /I.  We  then 
process  a  Sobel  filtering  on  the  scalar  part  to  produce  a  lu¬ 
minance  gradient.  Next,  we  apply  the  previous  mask  on  this 
gradient  to  keep  only  the  needed  luminance  information.  In 
fact  when  a  colour  vector  is  quite  collinear  to  the  greyscale 
axis,  the  saturation  gradient  given  by  the  previous  method 
is  not  appropriate.  Finally  this  is  merged  to  the  saturation 
gradient  in  order  to  also  detect  greyscale  variations.  We  see 
for  example  in  Figure  5  that  the  method  using  the  luminance 
information  merged  to  the  maximum  saturation  gradient  is 
more  efficient  as  there  are  more  details  on  the  left  parrot 
head.  There  are  as  well  big  achromatic  regions  behind  the 
parrots  detected  with  this  method  that  are  not  with  the  satu¬ 
ration  gradient  method  only. 

This  geometric  product  used  in  the  second  part  of  our  al¬ 
gorithm  could  have  been  processed  with  a  different  ju  vector, 
it  would  have  given  the  same  geometric  comparison  to  the 
colour  vector  chosen.  As  our  previous  gradient  was  based  on 
saturation,  we  needed  to  add  luminance  information  and  that 
is  why  we  chose  to  fill  /I  with  the  greyscale  axis. 

As  the  perception  of  colour  is  linked  to  the  Human  Visual 
System,  it  may  be  difficult  to  compare  results  in  colour  im¬ 
age  processing  with  greyscale  quantitative  evaluation  meth¬ 
ods.  We  could  have  compared  our  results  to  colour  images 
segmentation  databases  but  the  goal  of  this  paper  was  more 
to  show  how  using  Geometric  Algebra  could  help  to  perform 
a  colour  gradient.  We  can  nevertheless  compare  our  method 
visually  to  other  ones  in  Figure  6.  Di  Zenzo  method  [12]  is 
based  on  a  vectorial  approach  (c).  The  marginal  method  (b)  is 
a  fusion  of  three  different  gradients  done  on  the  three  colour 
channels  of  the  image.  It  is  quite  effective  on  this  picture. 
Results  show  that  our  method  detects  colour  edges  properly 
with  the  house  image  where  walls,  roof  and  sky  are  well  sep¬ 
arated  (e).  Because  the  luminance  in  (e)  is  considered,  we 
see  for  instance  the  details  of  the  window’s  house,  the  drain¬ 
pipe  and  the  bottom  of  the  roof  that  do  not  appear  on  the 
maximum  saturation  approach  (d).  The  Carron  method  (f) 
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Figure  5:  Colour  gradient  examples  :  (a)  original  image  ; 
(b)  colour  gradient  with  a  luminance  adjustment;  (c)  thresh¬ 
olding  on  the  colour  gradient ;  (d)  thresholding  on  the  maxi¬ 
mum  saturation  gradient  (without  the  modification  of  the  fil¬ 
ters). 


is  based  on  a  hue  marginal  gradient  but  it  is  completed  by 
luminance  and  saturation  gradients  when  hue  is  not  enough 
[1]. 

The  edges  obtained  by  our  approach  (d  and  e)  are  much 
thicker  than  in  all  the  other  methods.  The  Carron  approach 
gives  the  best  result  but  it  is  based  on  hue,  luminance  and 
saturation  variations  whereas  our  method  uses  only  the  two 
latter.  Our  method  gives  about  the  same  results  than  the  Di 
Zenzo  one  with  thicker  edges.  As  we  work  directly  in  the 
RGB  colour  space,  we  do  not  need  to  determine  a  geometric 
distance  between  colours;  we  use  a  saturation  distance  plus  a 
luminance  one  when  necessary. 

To  conclude,  we  propose  an  original  method  based  on 
geometric  algebra  to  detect  colour  edges.  This  formalism 
allows  to  manipulate  and  to  compare  the  colours  as  vectors 
in  the  3D  space.  This  method  gives  a  geometrical  gradient 
that  can  detect  colour  edges  properly  compared  to  already 
known  approaches  like  Di  Zenzo  or  Carron  ones.  It  uses 
RGB  colour  space  but  detects  saturation  and  luminance  vari¬ 
ations  with  geometrical  operations  allowed  by  ^3.  As  geo¬ 
metric  algebras  are  not  limited  to  three  dimensions,  an  exten¬ 
sion  to  N-dimension  from  this  work  is  possible.  Multiband 
images  as  those  produced  by  satellite  sensors  for  example  are 
being  studied  to  be  processed  by  geometrical  filters  in 

The  authors  thank  Poitou-Charentes  Region  board  which 
has  funded  this  work. 
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